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Abstract 

Background: Dugesio siculo is the only species of its genus not presenting an endemic or restricted distribution 
within the Mediterranean area. It mostly comprises fissiparous populations (asexual reproduction by body division 
and regeneration), most likely sexually sterile, and characterized by an extremely low genetic diversity interpreted as 
the consequence of a recent anthropic expansion. However, its fissiparous reproduction can result in an apparent 
lack of diversity within the species, since genetic variation within individuals can be as large as between them 
because most individuals within a population are clones. We have estimated haplotype and nucleotide diversity of 
cytochrome oxidase / within and among individuals along the species distribution of a broad sample of D. sicula, 
including asexual and the two only sexual populations known today; and predicted its potential distribution based 
on climatic variables. Our aim was to determine the centre of colonisation origin, whether the populations are 
recent, and whether the species is expanding. 

Results: The species presents 3 most frequent haplotypes, differing in a maximum of 1 1 base pairs. As expected 
from their fissiparous mode of reproduction, in half of all the analysed localities many individuals have multiple 
heteroplasmic haplotypes. The distribution of haplotypes is not geographically structured; however, the distribution 
of haplotypes and heteroplasmic populations shows higher diversity in the central Mediterranean region. The 
potential distribution predicted by climatic variables based modelling shows a preference for coastal areas and fits 
well with the observed data. 

Conclusions: The distribution and frequency of the most frequent haplotypes and the presence of heteroplasmic 
individuals allow us to gain an understanding of the recent history of the species, together with previous 
knowledge on its phylogenetic relationships and age: The species most probably originated in Africa and dispersed 
through the central Mediterranean. After one or multiple populations became triploid and fissiparous, the species 
colonized the Mediterranean basin, likely both by its own means and helped by human activities. Its present 
distribution practically fulfils its potential distribution as modelled with climatic variables. Its prevalence in coastal 
regions with higher water temperatures predicts a likely future expansion to northern and more interior areas 
following the increase in temperatures due to climate change. 

Keywords: Heteroplasmy, Fissiparity, Asexual, Recent colonization, Mediterranean basin, Trade routes, Niche 
modelling 



* Correspondence: mriutort@ub.edu 

Departament de Genetica, Facultat de Biologia and Institut de Recerca de la 
Biodiversitat (IRBio), Universitat de Barcelona, Barcelona, Catalonia, Spain 



O 



© 2013 Lazaro and Riutort; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the 
BiolVlGCl Central Creative Commons Attribution License (http://creativecommons.Org/licenses/by/2.0), which permits unrestricted use, 

distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public 
Domain Dedication waiver (http://creativecommons.Org/publicdomain/zero/1.0/) applies to the data made available in this 
article, unless otherwise stated. 



Lazaro and Riutort BMC Evolutionary Biology 2013, 13:268 
http://www.biomedcentral.com/1471 -21 48/1 3/268 



Page 2 of 14 



Background 

Dugesia sicula Lepori, 1948 is a freshwater planar ian 
(phylum Platyhelminthes, order Tricladida, suborder 
Continenticola) typically found in ponds, streams and 
springs close to the Mediterranean coast. The first rec- 
ord of D. sicula described a population comprising both 
sexual and asexual (fissiparous) individuals in Catania, 
Sicily [1]. Fissiparity is an asexual form of reproduction 
that involves transversal division of the individuals into 
two fragments and subsequent regeneration of the ab- 
sent structures. Fissiparous planarians do not develop a 
copulatory apparatus, which precludes the assignment of 
these populations to any species since this structure pro- 
vides most of the defining characters for species descrip- 
tion in planarians. However, the morphological characters 
of the copulatory apparatus of the sexual specimens 
allowed the species to be described. D. sicula was subse- 
quently found on the islands of Elba [2] and Mallorca [3]. 
None of the originally described sexual populations can be 
currently found [4]; but, some sexual populations in 
northern Africa (Tunisia and Algeria) have been recently 
reported [5-7]. Moreover, a molecular-based study [8] 
showed that most of the numerous fissiparous populations 
of genus Dugesia present in the western Mediterranean 
belong to D. sicula. Under certain conditions (mostly in 
the laboratory) triploid fissiparous planarians can develop 
hyperplasic ovaries and a copulatory apparatus [4,9,10], 
these individuals are known as ex-fissiparous. In D. sicula, 
it has been observed that some ex-fissiparous individuals 
laid cocoons but remained sterile [4,9-11], which indicated 
that asexual populations have lost the capability to sexu- 
ally reproduce. 

The literature suggests that D. sicula is widely distrib- 
uted in the Mediterranean area from Morocco to Greece 
[4,9,12], and our own observations show that this species 
mostly comprise fissiparous populations, most likely 
sexually sterile, distributed throughout the Mediterra- 
nean, as well as the Canary Islands, being the only spe- 
cies of this genus in the Mediterranean with a such a 
wide distribution. Contrasting with this wide distribution 
the species presents extremely low variability for the 
COI gene [8] among populations located hundreds of 
kilometres apart. These two features (wide distribution 
and low genetic variability) could be a consequence of a 
recent expansion promoted by human activities, as has 
been shown for other organisms [13]. Alternatively, the 
low levels of variability could be exclusively a conse- 
quence of their fissiparous reproduction, as described in 
Schmidtea mediterranea [14] and other organisms [15], 
and its wide distribution may be a consequence of its ac- 
tive spread over the years. In both cases, the hypotheses 
contrast with the poor dispersion capability assumed for 
freshwater planarians [16,17]. Freshwater planarians do 
not exhibit larval dispersal stages or forms resistant to 



desiccation; these individuals thus require contiguous 
freshwater bodies to survive and disperse [16,17]. How- 
ever, introductions due to human activities have been 
documented in other freshwater planarians, such as 
Girardia tigrina [18] and Schmidtea polychroa [19] and 
more recently have also been proposed for other Dugesia 
species [20], but in all these cases the introduced ani- 
mals either have a restricted distribution (Dugesia) or its 
introduction and progression in the new areas has been 
followed by scientist. 

Fissiparous reproduction can have effects on the distri- 
bution of genetic variability within the populations and 
even on individuals that could perhaps help disentangle 
this situation. Regeneration after fission is driven by neo- 
blasts [21,22], which are pluripotent stem cells. A new 
mutation in a neoblast would expand in the individual 
when this neoblast generates new tissues after the fission 
of the animal, thus increasing frequency of the mutation 
in the individual. Following replication via several fission 
cycles, the new mutation could disappear or spread in 
the population of cells within the individual and the or- 
ganism s population by genetic drift (Additional file 1: 
Figure SI). In this manner, we would expect to observe 
cells with different sequences within one individual. This 
could be considered a special case of heteroplasmy that 
differs from cases in which heteroplasmy is transmitted 
sexually or parthenogenetically through the mitochon- 
dria present in the oocyte (and the spermatozoon in 
some cases) where all cells from a new individual will 
present the same haplotypes. Heteroplasmy has been de- 
tected in many species [23], and some studies in humans 
have shown it increases with age in somatic cells and is 
more frequent in certain tissues, such as muscle tissue, 
most likely because of increases in the number of muta- 
tions caused by the presence of oxidative radicals 
[23,24]. Heteroplasmy has also been detected in COI se- 
quences of fissiparous populations of D. japonica [25]. 
Due to their characteristics, in fissiparous planarians the 
somatic heteroplasmy would be maintained in popula- 
tions and increased over time. 

If this situation holds in fissiparous populations of D. 
sicula, the analysis of their genetic variability could help 
in discriminating between a recent or ancient origin of 
the species populations and hence provide clues to the 
origin of its distribution. The absence or low frequency 
of heteroplasmic specimens in a specific locality of fis- 
siparous reproducing individuals could be interpreted as 
a consequence of a recent bottleneck, potentially be- 
cause of recent colonisation by a few individuals or re- 
covery from a recent population-wide decrease. The 
alternate condition, the presence of many heteroplasmic 
individuals, particularly those exhibiting a high number 
of different low- frequency haplotypes, could then poten- 
tially be related to aged individuals and populations, 
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although a recent colonisation by highly heteroplasmic 
individuals could also explain this outcome. Following 
the same reasoning, one would not expect to find het- 
eroplasmy, or different copies of mitochondrial genes 
within an individual in sexual organisms. 

In the present study, we analysed the genetic diversity 
and structure (using COI sequences) of a broad sam- 
pling of D. sicula from the Mediterranean, and the po- 
tential species' distribution was modelled to determine 
(1) whether intrapopulation and intraindividual nucleo- 
tide and haplotype diversity varies within the distribu- 
tion; (2) whether the species is currently expanding. 
With this data we expect to be able to answer questions 
on where the centre of origin of D. sicula s present dis- 
tribution is located, if the species has expanded recently 
assisted by humans and whether it is expected to con- 
tinue its expansion. 

Methods 

DNA extraction and sequence amplification 

Fifty-eight populations were sampled from the Mediter- 
ranean coast (Table 1 and Figure 1). Individuals from 
the Boussadia and Soukra populations (s39 and s36) in 
Tunisia reproduce sexually [5,6]. Individuals collected 
from Crete (s55), Greece, also possess a copulatory ap- 
paratus. Approximately 15 individuals (between 1 and 
20, Additional file 2: Table SI) were analysed from each 
location. High-molecular-weight DNA was purified from 
live or ethanol-fixed specimens using DNAzol reagent 
(Molecular Research Center Inc., Cincinnati, OH). A 
711-bp fragment of the mitochondrial gene COI was 
amplified and sequenced with specific primers (Tables 2 
and 3). In some individuals, the sequence had to be ob- 
tained in two non-overlapping fragments since we were 
not able to obtain amplifications in a single fragment. 
Sequencing was performed using Big Dye (3.1, Applied 
Biosystems, Norwalk, CT, USA) and an automated se- 
quencer, ABI Prism 3730 Applied Biosystems/Hitachi 
(Unitat de Genomica dels Serveis Cientifico-Tecnics de 
la UB) or at Macrogen Inc. (Korea). The DNA sequences 
were aligned in BioEdit [26] by sight based on the amino 
acid sequence and submitted to GenBank [GenBank: 
KC536630 to KC536644]. 

Intraindividual polymorphism 

To analyse the presence of multiple copies of the mito- 
chondrial genome within individuals, some of the PCR 
products were cloned using an HTP TOPO TA Cloning 
Kit for sequencing (Invitrogen, California, USA) follow- 
ing the manufacturers' protocol. We cloned the PCR 
products from sexual (5 individuals) and fissiparous indi- 
viduals, within the latter group we included both speci- 
mens who showed double peaks in the sequences 
obtained directly from the PCR product (11) and some 



that did not show this pattern (5). Approximately ten 
colonies from each individual were amplified and se- 
quenced using the T3 and T7 primers (included in the 
kit). Alignment was performed by sight based on the 
amino acid sequence, and the sequences were submitted 
to GenBank [GenBank: KC577271 to KC577351]. We 
tested whether the number of substitutions observed in 
the cloned sequences agreed with the expected random 
accumulation of mutations promoted by the polymerase 
errors to determine whether the variation observed in 
the sequences was caused by errors in the PCR process. 
Therefore, we assumed that errors in the PCR process 
followed a Poisson distribution, and the expected mean 
was computed using the polymerase error rate (6.3 x 10" 4 ) 
obtained by Savolainen et al. [28] and following their pro- 
cedure. We rejected the random accumulation of substitu- 
tions due to the occurrence of polymerase errors when 
the observed number of mutations fell outside the 99% 
confidence interval computed for the Poisson distribution 
(calculated with the DnaSP v5.10 program [29]). 

In a functional protein coding sequence, more changes 
are expected in the 1st and 3rd codon positions. To test 
the functionality of the sequences obtained in this study 
we calculated the expected number of changes per 
codon site in the case of a non-functional sequence 
(1:1:1) and tested whether the sequences obtained from 
the cloned PCR products fit the expectation, using a 
Chi-square test. 

Nucleotide and haplotype diversity analyses 

The DnaSP v5.10 program [29] was used to determine 
the haplotype and nucleotide diversity for each popula- 
tion and the species using only individuals that were not 
polymorphic, since we could not know the exact se- 
quence (phase) of haplotypes for all the individuals pre- 
senting multiple polymorphic sites. Those populations 
for which only one sequence was available were not in- 
cluded. Haplotype and nucleotide diversity were also 
calculated for each of the cloned individuals. Two haplo- 
type networks were constructed to study the relation- 
ships among haplotypes using TCS vl.21 [30]: (1) a 
haplotype network including only sequences from non- 
polymorphic individuals obtained after direct sequencing 
the PCR products and (2) a network containing only the 
cloned sequences from the polymorphic specimens. 

Distribution modelling 

To study the possible expansion of D. sicula species, 
their potential distribution was modelled with MAXENT 
software v3.3.3 k [31-33] using the default settings, a 
random test of 25%, the threshold rule of 10 percentile 
training presence, and 100 replicates. The independence of 
19 climatic variables were tested based on the R 2 statistic: an- 
nual mean temperature, mean diurnal range, isothermality, 
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Table 1 Sampling location details 



Code 


Origin 


Collectors and sampling date 


Latitude 


Longitude 


sOI 


PN Garajonay, La Gomera, Canary Is. 


E. 


Mateos, 2010 


28.1290 


-17.2515 


s02 


Setti Fatma, Ourika, Morocco 


A. 


H. Harrath, M. Yacoubi, 2009 


31.2192 


-7.6755 


s03 


Morocco 


P. Aguilera, C. Hernando, 1. Ribera, 2007 


33.5498 


-6.7486 


s04 


Andalusia, Spain 


M. Vila-Farre, 2006 


36.8070 


-5.3264 


s05 


River Chillar, Malaga, Spain 


N. Bonada, 2008 


36.7451 


-3.8769 


s06 


Es Corralassos, Eivissa, Balearic Is. 


C. Maritur, 2006 


38.9243 


1 .4362 


s07 


Algendar, Menorca, Balearic Is. 


S. 


Pons, 2006 


39.8730 


4.2090 


s08 


Binimel • la, Salairo, Menorca, Bal. Is. 


S. Pons, 2006 


40.0535 


4.0580 


s09 


Antella, Valencia, Spain 


M 


Vila-Farre, L. Leria, E. M. Lazaro, 2012 


39.0792 


-0.5926 


slO 


Casas del Rio, Valencia, Spain 


M 


Vila-Farre, L. Leria, E. M. Lazaro, 2012 


39.2967 


-1.1349 


si 1 


Jucar, Villalba de la Sierra, Cuenca, Spain 


M 


Vila-Farre, L. Leria, E. M. Lazaro, 2012 


40.2266 


-2.0894 


s12 


Ullals, Tarragona, Spain 


M 


Vila-Farre, L. Leria, E. M. Lazaro, 2012 


40.6912 


0.6960 


si 3 


Torre de la Carrova, Tarragona, Spain 


M 


Vila-Farre, L. Leria, E. M. Lazaro, 2012 


40.7527 


0.5666 


s14 


Duesaigues, Tarragona, Spain 


M 


Vila-Farre, L. Leria, E. M. Lazaro, 2012 


41.1516 


0.9283 


si 5 


Tres Pins Nursery, MontjuTc, Barcelona, Spain 


M 


Riutort, E. M. Lazaro, 2006 


41.3676 


2.1620 


s16 


MontjuTc, Barcelona, Spain 


M 


Vila-Farre, 2012 


41.3677 


2.1622 


si 7 


P. Laberint, Barcelona, Spain 


M 


Vila-Farre, 201 1 


41.4400 


2.1451 


s18* 


Orfes, Girona, Spain 


M 


Riutort, 1996 


42.1717 


2.8679 


s19 


Fluvia, Girona, Spain 


M 


Riutort, 201 1 


42.1610 


2.9585 


s20 


Le Pont de Reynes, Languedoc-Roussillon, France 


M 


Vila-Farre, L. Leria, E. Sola, 201 1 


42.4958 


2.7148 


s21 


Lunac, Languedoc-Roussillon, France 


M 


Vila-Farre, L. Leria, E. Sola, 201 1 


43.7082 


3.1955 


s22 


Prades-le-Lez, Languedoc-Roussillon, France 


M 


Vila-Farre, L. Leria, E. Sola, 201 1 


43.6841 


3.8605 


s23 


River Herault, France 


M 


Riutort, E. M. Lazaro, 2009 


43.9080 


3.7371 


s24 


Invrea, Liguria, Italy 


M 


Pala, 2011 


44.3831 


8.6096 


s25 1 


Bisagno, Liguria, Italy 


M 


Pala, 2011 


44.4477 


9.001 1 


s26 


Rapallo, Genova, Italy 


M 


Pala, 2011 


44.3522 


9.2308 


s27* 


Isola de I'Asinara, Sardinia, Italy 


M 


Pala, 2005 


41.0707 


8.3198 


s28* 


Isola de Tavolara, Sardinia, Italy 


M 


Pala, 2005 


40.9093 


9.6093 


s29* 


Torre Argentina, Sardinia, Italy 


M 


Pala, 2005 


40.3282 


8.4610 


s30* 


Acqua sa Murta, Sardinia, Italy 


M 


Pala, 2005 


39.0768 


8.4238 


s31* 


Acqua sa Canna, Sardinia, Italy 


M 


Pala, 2005 


39.0638 


8.4543 


s32 


River Joumine, Tunisia 


R. Sluys, A. H. Harrath, 2006 


36.9871 


9.6038 


s33* 


Fourna, Tunisia 


M. Charni, 2007 


36.6935 


9.8306 


s34* 


Lebna, Tunisia 


M. Charni, 2007 


36.7385 


10.9220 


s35* 


Kherba, Tunisia 


M. Charni, 2007 


36.8822 


10.9041 


s36 2 


Soukra, Tunisia 


M. Charni, 2007 


36.8685 


10.2656 


s37* 


r~hih^ Ti inKi^ 
1 1 ua, i u 1 1 loia 


M. Charni, 2007 


36.6970 


1 0.771 8 


s38* 


Siliana, Tunisia 


M. Charni, 2007 


36.0822 


9.3878 


s39 2 


Boussadia, Tunisia 


M. Charni, 2007 


36.1011 


9.6118 


s40* 


Lakhmes, Tunisia 


M. Charni, 2007 


35.9978 


9.4813 


s41 


ATn Ghanim, Tunisia 


A. 


H. Harrath 


36.0589 


10.4025 


s42 


Leben, Tunisia 


A. 


H. Harrath, 2007 


34.5801 


9.8849 


s43 3 


Gabes, Tunisia 


R. Sluys, A. H. Harrath, 2006 


34.1739 


10.0133 


s44 


ATn Chebika, Tunisia 


M. Riutort, 2008 


33.9283 


8.1235 
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Table 1 Sampling location details (Continued) 



S4j 


Fra. Di Marsala, 0 Sossio, Sicily 


M. Pala, G. A. Stocchino, E. 


M. Lazaro, 2009 


37 77QA 
J/.//O0 


11/1 ^Q^ 


s46 


Marausa, Sicily 


M. Pala, G. A. Stocchino, E. M.Lazaro, 2009 


Q 7 fio^n 

o/.yooy 


1 2.5345 


54/ 


Sorgente, Sicily 


M. Pala, G. A. Stocchino, E. 


M. Lazaro, 2009 




I z.yuiQ 


r/l Q 

S4o 


Castellammare del Golfo, Sicily 


M. Pala, G. A. Stocchino, E. 


M. Lazaro, 2009 


OO.Uz I D 


i z.yuoy 


s49 


Ponte Saraceni on river Simeto, Sicily 


/VI . raid, ZUUD 




37.701 2 


1 a onn 1 
1 4.oUU I 


s50 


Casale del Simeto, Sicily 


M. Pala, 2007 




37.5804 


1 4.8737 


s51 


Smertos-Filiates, Greece 


E. Sola, K. Gritzalis, 2010 




39.6196 


20.2537 


s52 


Cephalonia, Greece 


R. Sluys, 2009 




38.1591 


20.7263 


s53 


Kirinthos-Mandoudi, Euboea, Greece 


E. Sola, 2010 




38.8024 


23.4473 


s54 


Nomia-Kalives, Peloponnese, Greece 


E. Sola, J. Sola, 2010 




36.6555 


22.9956 


S55 4 


Crete, Greece 


E. Sola, E. Mateos, 2009 




35.0803 


25.6508 


s56 5 


Rhodes, Greece 


E. Sola, E. Mateos, 2009 




36.3371 


28.0625 


s57 


Samos, Greece 


M. Vila-Farre, 2010 




37.7767 


26.7258 


s58 6 


Chios, Greece 


M. Vila-Farre, 2010 




38.5273 


25.8705 



All populations were identified as D. sicula by their COI sequence except populations marked with * that were previously identified by karyotype; 1 Cohabited with 
D. liguriensis; 2 Sexual populations with karyotype 2n = 18; 3 Collected in the same pond where D. maghrebiana has been reported; 4 Individuals with copulatory 
apparatus; 5 Cohabited with D. elegans; 6 Exfissiparous. 



temperature seasonality, maximum temperature of 
warmest month, minimum temperature of coldest month, 
temperature annual range, mean temperature of wettest 
quarter, mean temperature of driest quarter, mean 
temperature of warmest quarter, mean temperature of 
coldest quarter, annual precipitation, precipitation of wet- 
test month, precipitation of driest month, precipitation 
seasonality, precipitation of wettest quarter, precipitation 
of driest quarter, precipitation of warmest quarter and 
precipitation of coldest quarter. Finally four independ- 
ent climatic variables were used: isothermality, mean 
temperature of the wettest quarter, mean temperature 
of the driest quarter and precipitation seasonality. From 



61 known localities (58 analysed in this study and 3 se- 
lected from the literature [7] -Algeria- and personal 
communications of L. Leria -Mallorca- and E. Sola - 
North of Italy-), 48 observation points were used, and 
13 locations were discarded because of their geographic 
proximity (populations situated within the same cell of 
the climate layer: sl3, sl5, sl7, sl9, s20, s22, s25, s30, 
s32, s34, s35, s45 and s49). 

Results 

Population analyses of non-polymorphic individuals 

In total, for the 58 locations sampled (Table 1), 619 se- 
quences were obtained, and 163 (Additional file 2: Table SI) 





slSOO 




sl9 




sl7 


sl4 

O 

sl2 

QOsl3 


sl6(9 sl5 



sllQ 



o 

s05 




O s03 




o 

s44 



9 s ' 



9 s49 




O 


J <£)s50 




s52 


s45 






s32 


s36 


s35 


• 


O 




s33# 




7 s34 




S3" 


^o6 39 0 






s41 


s40 






Specimens 


from 


s39 



& i 

# s56 

s55 



to 



reproduction. The individuals from s55 possessed a copulatory apparatus. s43 was collected in the same pond where D. maghrebiana has 
been reported. s56 cohabited with D. elegans, and s25 cohabited with D. liguriensis. Details regarding the sampling sites and collectors are 
provided in Table 1 . 
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Table 2 Primers used for amplification and sequencing 

Name Sequence 5-3' An. Temp. fC) 1 Source 

BarT (F)* ATGACDGCSCATGGmAATAATGAT 43 Alvarez-Presas et ol., [27] 

COIF (F) CCNGGDmGGDATDR^A^CWCA 45 Lazaro et ol. [8] 

COIR (R) CCWGTYARMCCHCCWAYAGTAAA 45 Lazaro et ol. [8] 

COIbarc_plat_R (R) TAATOAAATATAAACCTCAGGATG 40 Lazaro et ol. [8] 

BARSI (F)** GATAGGCGG^GGTAAATGG 50 This study 

F forward, R reverse, *Only for amplification, **Only for sequencing 1 Annealing Temperature. 



of these sequences showed double peaks in the chromato- 
gram and were discarded for the population analyses. Be- 
cause some populations had only one sequenced specimen 
(s04, si 8, and s30) or all of their individuals exhibited poly- 
morphic positions, only 453 individuals from 51 populations 
could be analysed. Only 16 populations showed more than 
one haplotype, and the majority of these showed two 
haplotypes. 

The haplotype network was based on an alignment of 
456 sequences and a length of 624 base pairs. The total 
number of haplotypes for the COI gene in the popula- 
tions was 15 (Figure 2, Additional file 2: Table SI). The 
most frequent haplotype (haplotype A, the red circle in 
Figure 2) was observed in 292 individuals from 41 popu- 
lations throughout the distribution, and presented 11 
differences regarding haplotype B (the second most fre- 
quent, blue circle in Figure 2). Haplotype B was present 
in 144 individuals from 18 populations, the majority of 
which were found northwest of the Mediterranean. 
Haplotype C (brown in Figure 2) was only observed in 
six individuals from southwestern populations, and pre- 
sented six differences regarding haplotype A. The other 12 
haplotypes were unique or observed in only two individ- 
uals. Haplotypes A and B were observed simultaneously in 
6 populations from the centre of the Mediterranean. 

Intra-individual polymorphism 

Between one and ten polymorphic positions (double 
peaks in the sequences) were observed in 163 individuals 
from 29 populations (one-half of the populations, the 
black circles in Figure 3) when sequencing the PCR 
product directly. Twenty-one PCR products (11 from in- 
dividuals exhibiting sequences with double peaks and 10 
without) from 13 populations were cloned, and ap- 
proximately 10 clones per individual (Table 4 and 
Additional file 2: Table S2) were sequenced to determine 



Table 3 Length of amplified fragments 



Primers 




Sequence length 


BarTVBARSI** - 


hCOIR 


711 


BarTVBARSI** - 


hCOIbarc_plat_R 


396 


COIF + COIR 




228 



*Only for amplification, **Only for sequencing. 



the haplotypes of those individuals. The final alignment 
had 200 sequences and a length of 537 bp. 

In all cases we found multiple haplotypes within each 
individual even for those not showing polymorphic sites 
in the direct PCR sequence. However, some of the point 
mutations detected in the cloned sequences could have 
been an artefact of the PCR and/or cloning processes. 
Savolainen et al. [28] estimated this error rate to be 6.3 x 
10" 4 . We compared the number of expected point muta- 
tions based on this rate with the number observed in each 
cloned specimen, taking into account the entire length of 
the amplified and sequenced DNA (the number of clones 
x the number of positions; Table 4). All of the individuals 
with polymorphic sites exhibited significant differences 
between the two values (except for the s48.10), which indi- 
cated that the detected point mutations could not be 
accounted for by the addition of erroneous nucleotides by 
the polymerases during the PCR reactions. None of the 
sexual individuals analysed had a significant result. More- 
over, the cloning results showed that haplotypes A, B and 
C were present simultaneously in some of the specimens 
(Additional file 2: Table S2). However, three of the speci- 
mens in which we had not detected polymorphic sites 
when sequencing directly from PCR, also showed signifi- 
cant differences between the expected and observed point 
mutations (Table 4). This result suggests that, in these in- 
dividuals, despite that some point mutations could be 
caused by polymerase error, there are most likely low- 
frequency haplotypes that we cannot detect in the direct 
sequences using PCR. In the subsequent analyses, we con- 
sidered all of the individuals that showed a significant re- 
sult to be heteroplasmic and the remaining individuals to 
be non-heteroplasmic. 

For the haplotype network reconstruction (Figure 4 
and Additional file 1: Figure S2), the single haplotypes 
sequenced from the non-heteroplasmic individuals were 
discarded since those must most probably be the result 
of a polymerase error as shown in the previous tests. 
The final network contained 129 sequences and 57 dif- 
ferent haplotypes were observed; haplotypes A, B and C 
were the most frequent and appeared 46, 43 and 16 
times, respectively. The remaining haplotypes were 
unique in the majority of the cases and showed a single 
difference from the A, B or C haplotypes. All of the 



Lazaro and Riutort BMC Evolutionary Biology 2013, 13:268 
http://www.biomedcentral.com/1471 -21 48/1 3/268 



Page 7 of 14 



4)S G 





Figure 2 Haplotype network and distribution in the studied populations. Red corresponds to haplotype A, dark blue corresponds to 
haplotype B, and brown corresponds to haplotype C. Only haplotypes from non-polymorphic individuals are represented. The grey dots on the 
map show populations in which all of the specimens showed polymorphic sites. 



analysed individuals, except specimen s33.10 (Figure 4 
and Additional file 2: Table S2), exhibited at least one 
of the most frequent haplotypes (A or B), and all of the 
individuals showed at least one unique haplotype 
(Additional file 1: Figure S2). The haplotype and nu- 
cleotide diversity, number of polymorphic sites and 
number of haplotypes were higher in heteroplasmic 
than non-heteroplasmic populations (H D = 0.847 vs. 
0.624, tt = 0.0096 vs. 0.0017, S= 14.1 vs. 3.9 and h = 6.2 
vs. 4.5, Additional file 2: Table S2). 

The values of n for heteroplasmic individuals (0.0096) 
is slightly higher than that for all populations consider- 
ing only non-polymorphic individuals (0.0081), and in 



general intra- individual n values are as high or even 
more than the n value found for the whole species. For 
H Dj there is an increase when we only consider the het- 
eroplasmic individuals due to the presence of multiple 
haplotypes at similar frequencies. 

Predicted distribution 

The MAXENT results show the potential distribution of 
D. sicula from the Mediterranean coasts, extending 
throughout the western region: northern Tunisia, 
Algeria and Morocco, southern Iberian Peninsula, the 
Italic Peninsula, and islands of the Balearic archipelago, 
Sicily and Sardinia. All of the populations were observed 




Figure 3 Locations with polymorphic individuals. The black dots represent locations with polymorphic individuals. The white dots represent 
populations in which all of the individuals are non-polymorphic. 
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Table 4 Polymerase mistake versus heteroplasmy test 



Specimen 


Direct sequence 
(heteroplasmic positions) 


N cloned 
sequences 


N bp 


Expected point 
mutations 


Observed point 
mutations 


Poisson 
P[X < obs] 


S23.02 


Haplotype A 


8 


4296 


2.7 


1 


0.07 


S23.09 


Haplotype A 


7 


3759 


2.3 


1 


0.10 


S27.14 


Haplotype B 


11 


5907 


3.7 


14 


1.00* 


S27.19 


Haplotype B 


11 


5907 


3.7 


11 


1.00* 


S36.12 


Haplotype B 


10 


5370 


3.4 


7 


0.95 


S36.03 1 


Haplotype A 


10 


5188 


3.0 


4 


0.64 


S39.08 


Haplotype A 


7 


3759 


2.3 


6 


0.97 


s39.11 


Haplotype A 


10 


5370 


3.4 


7 


0.95 


S39.14 


Haplotype B 


10 


5370 


3.4 


4 


0.57 


S42.10 


Haplotype B 


9 


4833 


3.0 


34 


1.00* 


S06.07 


Heteroplasmic (2) 


8 


4296 


2.7 


33 


1.00* 


S26.12 


Heteroplasmic (1) 


10 


5370 


3.4 


31 


1.00* 


S31.18 


Heteroplasmic (10) 


29 


15530 


9.7 


219 


1.00* 


S31.20 


Heteroplasmic (6) 


9 


4833 


3.0 


63 


1.00* 


S33.10 


Heteroplasmic (6) 


4 


2148 


1.3 


5 


0.99* 


S47.03 1 


Heteroplasmic (3) 


9 


4692 


2.7 


71 


1.00* 


s47.11 1 


Heteroplasmic (5) 


10 


5088 


2.7 


32 


1.00* 


S48.04 


Heteroplasmic (4) 


12 


6444 


4.0 


95 


1.00* 


S48.10 


Heteroplasmic (2) 


10 


5370 


3.4 


4 


0.57 


S55.10 


Heteroplasmic (4) 


3 


1611 


1.0 


8 


1.00* 


S56.15 


Heteroplasmic (7) 


7 


3759 


2.3 


26 


1.00* 



Expected number of point mutations for each individual based on a polymerase error rate of 6.3 x 10~ 4 [28], point mutations found in cloned sequences, and 
probability that observed point mutations exceed the expected based on a Poisson distribution. 

Sequence length was 537 bp, except in some cases marked with \ N bp, total length amplified (length amplified x number of clones); ^Significant differences 
based on 99% criterion. 



within the predicted area, except for population s44 
(south of Tunisia) and sOl (Canary Islands), (Figure 5). 

Discussion 

Polymorphic sequences within individuals 

Heteroplasmy, the presence of different copies of the 
mitochondrial genome within the cells of an individual, 
has been reported for several taxa [24,34-39] and was an 
expected outcome in the case of fissiparous flatworms, 
where somatic heteroplasmy would be inherited through 
clonal reproduction (Additional file 1: Figure SI). How- 
ever, there are multiple hypotheses other than hetero- 
plasmy that can explain polymorphic sequences; these 
explanations include the presence of numts or errors in- 
troduced by polymerases during the PCR reaction. In 
the present case, the statistical analysis (Table 4) indi- 
cated that the number of mutations observed in some 
individuals cannot be explained by errors introduced by 
the polymerases during PCR. Moreover, although the full 
evidence for a mitochondrial sequence to be functional 
and not present in the nucleus (numts) can be reached 
only by sequencing mRNAs (cDNAs) from the individuals 



[36], in our case a variety of evidence makes us think 
that the diverse sequences found within single planar- 
ians are not numts. Numts typically have short se- 
quences (<200 bp) and their mutation percentages are 
between 5 and 25% [40-43]. The length of the shortest 
amplified fragments in this study was 228 bp, and the 
polymorphic positions observed in these fragments 
were also present when we amplified the larger frag- 
ment (711 bp) (Table 3). The percentage of variation 
for our polymorphic sequences was below 2.5%. Most 
of the substitutions were synonymous, and no stop co- 
dons were observed. Moreover, it is expected that the 
distribution of substitutions for a functional sequence 
shows a bias with more differences observed at the 
third and first position than at the second. We tested 
the distribution of changes in the cloned sequences of 
both those individuals in which the previous test could 
not reject that the substitutions seen are due to Taq 
error (non-heteroplasmic) and in those where the en- 
zyme cannot explain the number of changes observed 
(heteroplasmic). The result shows that in the latter 
there is a significant deviation from the 1:1:1 expected 
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Figure 4 Cloned specimen haplotype network and distribution. The pie charts represent the haplotypes observed in each cloned 

heteroplasmic individual (Additional file 2: Table S2). The cloned individuals that were not heteroplasmic do not have pie charts. Red corresponds 

to haplotype A, dark blue corresponds to haplotype B, and brown corresponds to haplotype C Colours of populations are as in Figure 1. 
\ ' ) 



ratio (Additional file 2: Table S3). However, some cases 
of recent numts have been described [36] that still pre- 
sented the characteristics of a functional gene (were 
not short or did not present high levels of variation). 
There is still one more thing in our case that provides 
support for all the copies being functional, in the 



cloned individuals presenting heteroplasmy, we found 
copies of at least two or the three most frequent haplo- 
types (A, B and C) that are also found alone in some in- 
dividuals and in whole populations, and that do not 
present stop codons, no traces of being non-functional 
(Additional file 2: Table S2). 




Figure 5 Potential distribution of D. sicula in the Mediterranean Basin. The green colour scale shows the predicted probability that 
conditions are suitable for the species in the studied area under the 10 percentile training presence criterion. The red dots show the observation 
points used for the prediction analysis, and the yellow dots show known localities discarded from the analysis (si 3, si 5, si 7, si 9, s20, s22, s25, 
s30, s32, s34, s35, s45 and s49). Localities surrounded by purple circle were found outside of the predicted area. 

V J 
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Although fission could easily explain the coexistence 
of low-frequency haplotypes that differ in one or a few 
substitutions in a given specimen (Additional file 1: 
Figure SI), the origin of heteroplasmic individuals that 
simultaneously present haplotypes A, B and C is diffi- 
cult to envision. They can be a consequence of past 
sexual reproduction episodes, or a result of the appear- 
ance of the third haplotype in individuals that already 
presented the other two. In any case, the finding of 
more unique haplotypes (Figure 4 vs Figure 2), differing 
mostly by one nucleotide from haplotypes A, B or C, in 
the heteroplasmic individuals supports the prediction 
that asexuals will tend to accumulate rare variants in 
all loci, which will be more accentuated in older asex- 
uals ([44] and references there in). 

On the other hand, as expected, heteroplasmy is not 
observed in sexual individuals. Although sexual individ- 
uals can also reproduce by fission, gametogenesis and 
zygote formation could represent a limitation for hetero- 
plasmy propagation, which is a situation also observed 
in plants that reproduce both sexually and asexually 
[23,35]. However, some specimens from the s55 popula- 
tion from Crete possessed a copulatory apparatus and 
were consequently considered sexual individuals. The 
observation that 3 out of the 13 sequenced individuals 
from this population were heteroplasmic could indicate 
that the individuals with a copulatory apparatus were ac- 
tually ex-fissiparous or that the population is composed 
of a mix of sexual and fissiparous individuals. 

Origin and expansion 

A recent study [8] emphasised that D. sicula distribution 
could be a consequence of a relatively recent dispersion 
driven by human activities in accordance with a previous 
report [45]. In the present study we found that in D. 
sicula populations, only two haplotypes showed high fre- 
quencies, differing at 11 sites over 624 bp, and most of 
the populations had a single haplotype (Figure 2 and 
Additional file 2: Table SI). Differing from what was 
found in S. mediterranea [14], a non- overlapping geo- 
graphical distribution for the different haplotypes, the 
most frequent D. sicula haplotypes are found practically 
throughout the distribution of the species. This lack of 
structure seems to reinforce the idea of a recent expan- 
sion. Moreover, the haplotype and nucleotide diversity 
values are extremely low at the intra-population level 
and when the COI values for the entire set of popula- 
tions were compared to those observed in other planar- 
ian species (Table 5 and Additional file 2: Table SI), such 
as the freshwater Schmidtea mediterranea [14] or the ter- 
restrial Cephaloflexa bergi and Imbira marcusi [27] . How- 
ever, the presence of numerous heteroplasmic individuals 
bearing rare haplotypes in many fissiparous populations 
found in the present study (the intraindividual variability 



Table 5 Comparison of diversity parameters between 
different triclad species 



Species 


H D 


71 


References 




Dugesia sicula 


0.491 


0.0081 


This study 




Schmidteo mediterranea 


0.773 


0.0223 


Lazaro et a I., [14] 




Cephaloflexa bergi 


0.950 


0.0573 


Alvarez-Presas et al. 


[27] 


Imbira marcusi a 


0.857 


0.0463 


Alvarez-Presas et al. 


[27] 



a this species has recently been renamed (Carbayo et al. [46]), in Alvarez-Presas 
et al. [27] appears as Geoplana goetschi Marcus, 1951. 



can be as high as the total variability of the species) 
suggests the ancient age of these lineages and likely for 
the populations where we find multiple heteroplasmic 
specimens. 

Similar to other long-lived asexual lineages, D. sicula 
follows the expected accumulation of neutral or slightly 
deleterious substitutions over time and independent of 
one another, which results in high heterozygosity in their 
populations [47,48]. However, in the case of D. sicula 
this accumulation is found within individuals, which 
does not allow us to exactly quantify the levels of intra- 
population variability (unless we cloned and sequenced 
many individuals from each population). However, tak- 
ing into account our predictions of the expected levels 
of heteroplasmy and of heteroplasmic individuals in 
populations as related to their age, we can draw some 
conclusions on the origin and distribution of the species. 

The high frequency and allocation of haplotypes A and 
B appears to indicate their antiquity in the species, par- 
ticularly haplotype A, which is observed from Morocco to 
Greece (Figure 2). In the northwestern populations 
(northern Catalonia and France), all 97 sequenced individ- 
uals showed haplotype B (one individual was polymorphic 
but likely non-heteroplasmic). The low genetic variability 
in this region could indicate a relatively recent colonisa- 
tion of this area by individuals deriving from geographic- 
ally distant populations situated in the central or eastern 
regions (no close population shows haplotype B, Figure 2). 

Haplotype C was observed only in western populations 
and is most likely derived from haplotype A based on 
the haplotype network (Figures 2 and 4, and Additional 
file 1: Figure S2). Haplotype C is observed in non- 
heteroplasmic individuals in Morocco and the southern 
Iberian Peninsula and may have spread from there to the 
Canary Islands, where C was the only haplotype ob- 
served. Haplotype C was not present in any sexual speci- 
men, but it was observed in heteroplasmic individuals 
from the western Mediterranean islands and northern 
Italy (Figure 4). Based on these data, we cannot deter- 
mine with certainty where the haplotype originated 
from, but we can assume it had a relatively more recent 
origin than the other two haplotypes (A and B) and dis- 
persed only through the western region. However, we 
can speculate that the presence of haplotype C in 
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heteroplasmic individuals in the central regions could 
indicate its origin is there, and its prevalence in the sOl, 
s02 and s05 populations could be the result of a founder 
effect during colonisation. 

Insufficient information provided by the few most 
common haplotype distributions and the fact that the 
presence of multiple single haplotypes in a population 
not necessarily are a consequence of an old age for the 
population at that locality, make it difficult to determine 
where colonisation began in the Mediterranean. How- 
ever, some facts can help: (1) the evolutionary sister 
group of this species, D. aethiopica [8] is distributed ex- 
clusively in Africa [49], and their split is dated around 4 
My (c. 7-2 My) [20], well after the Mediterranean basin 
was formed and all islands had reached their present lo- 
cation [50,51], (2) most of the populations (either sexual 
or fissiparous) in which haplotypes A and B cohabit are in 
the central region of the Mediterranean basin (Figure 2), 
(3) known sexual populations originate only in either Si- 
cily (now lost) or Tunisia and Algeria. Collectively, these 
data could indicate that the central Mediterranean is the 
region from where dispersal took place. However, the 
question as to how this species attained such a wide distri- 
bution with such a low dispersal capability [16,17] re- 
mains. Various factors could have contributed, including 
their fissiparous reproduction and their capability to resist 
high water temperatures. 

It has been shown [52] that asexual reproduction in 
planarians allows a more rapid increase in population 
size under ecological conditions in which food is limited. 
Moreover, a study by Vila-Farre [53] have shown that D. 
sicula is sensitive to low temperatures and to abrupt 
temperature changes, while Charni et al. [10] found that 
the diverse Tunisian populations lived in water with 
temperatures ranging from 13 to 22°C. It has also been 
observed that during spring and autumn, higher rates of 
fissioning occur [54]. Therefore, they can resist high 
temperatures of water, and can reproduce and increase 
their populations very fast. 

But could the species by itself reach its present distri- 
bution area? Given the maximum age of the species situ- 
ated around 4 My (c. 7-2 My) [20], when the 
Mediterranean basin was already formed and its islands 
in its present situation, it seems difficult to explain how 
the animals could have moved between islands and the 
continent (Corsica, Sicily, Tunisia...) even in one direc- 
tion or the other. If the origin of the species preceded 
the Messinian Salinity Crisis or else as a consequence of 
the sea level changes during the Pleistocene glaciations, 
we could maybe hypothesize the contact among some 
fluvial basins in the central Mediterranean that would 
allow its expansion there, but still this would not explain 
its presence in places situated as far away as the Greek 
islands, the Canary islands or how haplotype B could 



have reached the French coast. These latter cases pro- 
vide support for a mixed history of self-dispersal and hu- 
man introductions. De Vries [55] already proposed that 
given the antiquity of human trade in the Mediterranean 
basin there may have been multiple opportunities for 
these animals to be transported around, including mod- 
ern commercial trade of fish, plants or timber, as has 
been proposed for other planarian species [56-58], which 
also could explain this expansion. 

Our hypothesis is then that D. sicula most likely origi- 
nated in northern Africa and dispersed through the 
Central Mediterranean region. After one or multiple 
populations became triploid and fissiparous, they colo- 
nized the whole Mediterranean basin either by their own 
means or assisted by human activities (probably both 
scenarios occurred). 

Species distribution and future expansion 

Under the hypothesis that D. sicula is a relatively recent 
coloniser, it was advisable to test whether this species 
had reached its current range of fitted distribution or 
was still expanding. The potential distribution predicted 
with the climatic models for D. sicula (Figure 5), fits well 
with what has been observed, mainly occupying coastal 
regions. It has been reported that D. sicula does not 
cope well with abrupt changes in temperature [53] and 
is typically observed near the coast, where the water 
temperature is highest in the winter and the thermal 
changes are less abrupt than in the high mountains. All 
known populations, except for s44 and sOl, have been 
observed within the predicted area; however, population 
s44 inhabited an oasis in the Tunisian desert, the special 
conditions in this reduced area could be ideal for the 
species although precision of climate layers used in the 
analysis do not allow for detecting them. Population sOl 
in the Canary islands is found within Garajonay National 
Park, in a area of laurel forest, thought to be the remnant 
of subtropical woods that covered the Mediterranean area 
during the Tertiary period, which likely provides D. sicula 
the conditions needed to survive, although the models, 
again as a consequence of the scale, may have not been 
able to detect it. 

Some regions that are potentially suitable for D. sicula 
have not yet been widely sampled, such as northern 
Africa (except for Tunisia and some regions of Algeria 
and Morocco), the southwestern Iberian Peninsula and 
parts of the Italic Peninsula. However, restriction pattern 
analyses showed similarities among some populations 
north of the Italic Peninsula and D. sicula [59,60], al- 
though the Italian populations could not be identified; 
moreover, Benazzi and Deri [61] reported populations in 
Tuscany (including the island of Pianos), Ponza Island 
(between Rome and Naples), Rome (Lazio) and Calabria 
that could be attributed to D. sicula based on 
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morphological features, and some populations with a 
chromosomal complement of 27+ 2B chromosomes have 
been reported in Calabria by Deri et al [62]. In other 
areas, D. sicula has never been reported (Corsica) or is 
relegated to the periphery (Sardinia), although these re- 
gions appear to fulfil the climatic requirements of the 
species. However, several populations of D. benazzii have 
been observed in these regions. The competition be- 
tween D. sicula and D. benazzii may explain why D. 
sicula is absent in Corsica and has a low presence in 
Sardinia. However, D. sicula is present together with 
other planarian species in localities such as Lebna in 
Tunisia {S. mediterranea), which points to the possibility 
that historical reasons could also explain its absence in 
some regions (they may have never been introduced by 
humans or had the opportunity to reach by their own 
means). 

For the eastern Mediterranean, some populations 
described in Israel as Type 4 by Bromley [63] showed 
karyotypes and morphological characteristics similar to 
those of D. sicula. Although Type 4 was related to D. 
biblica in a study by Bromley [64], these two species 
cannot be definitively differentiated [12]. An accurate 
study of these two species is necessary to clarify their 
status and confirm whether D. sicula is also present in 
that region or whether the presence of a closely related 
species potentially precluded its colonisation. 

Recent exhaustive samplings performed in the Iberian 
Peninsula in search of Dugesia subtentaculata (M. 
Vila-Farre and L. Leria personal communication and 
one of the authors, MR), allows us to confidently state 
that D. sicula is only observed on the east and south 
coast, as the prediction analysis suggests. In this area, 
D. sicula cohabits with other freshwater planarian spe- 
cies such as S. mediterranea. In recent years, popula- 
tions of S. mediterranea [14] appeared to have 
vanished, whereas populations of D. sicula became 
more frequent. Additionally, in the Greek population 
from Rhodes (s56), D. sicula cohabits with the autoch- 
thonous D. elegans, and only a few animals belong to 
this latter species. This result could indicate that D. 
sicula can have an adverse effect on the native triclad 
species, which contrasts with what appears to occur in 
Sardinia. A similar situation has been reported in 
Wales, where the presence of two immigrant species, 
Planaria torva and Girardia tigrina has been followed 
via three surveys over 50 years [56,58,65,66]. The most 
recent report [57] indicates that the invasion of G. 
tigrina has been more successful, which is most likely 
because of its asexual reproduction and opportunistic 
feeding characteristics, although other factors, such as 
its capability to adhere to surfaces better than other tri- 
clads, have also been considered. Moreover, G. tigrina 
has been shown to have a negative effect on native 



species of triclads, which is most likely explained by 
food competition. 

Conclusion 

In conclusion, it appears that D. sicula has reached a large 
proportion of the area of its potentially favourable distri- 
bution in the Mediterranean basin, being a remarkable 
case of a broad colonisation, in extreme contrast with the 
rest of Mediterranean Dugesia species, with all of them 
being endemic or with very restricted distributions. D. 
sicula expansion is now limited to spreading to new fresh- 
water basins within the areas it currently inhabits. How- 
ever, future changes increasing the temperature, such as 
those predicted by climate change hypothesis, could ex- 
pand its fitted area to more northern and interior areas. 

Additional files 



Additional file 1: Figure SI. Schematic representation of the spreading 
of a new mutation in fissiparous organisms. The yellow and orange colours 
represent the tissues bearing the new mutations. The arrows represent fission 
(black) and regeneration (blue) cycles. Figure S2. The haplotype network 
based on cloned sequences. Each circle represents a haplotype observed in 
the sample. The black dots represent intermediate (non-present) haplotypes, 
the lines connecting the haplotypes (either present or not) represent one 
nucleotide change, and the size of each circle is proportional to the 
haplotype frequency in the sample. The colours in the network correspond 
to those on the map and represent the studied locations. The crossed-out 
names of the individuals represent non-heteroplasmic specimens, and their 
sequences were not included in the network. 

Additional file 2: Table SI. Genetic diversity in the populations. 

Table S2. Genetic diversity in the cloned individuals. Table S3. Chi-square 
test for the distribution of substitutions along the three codon positions 
(null hypothesis equifrequent distribution, 1:1:1). 



Competing interests 

The authors declare that they have no competing interests. 
Authors' contributions 

MR and EML designed the study. EML performed the molecular work and 
obtained the sequence data. EML and MR conducted all of the analyses and 
wrote the manuscript. Both authors read and approved the final manuscript. 

Acknowledgements 

This research was supported by DGI-MEC Grant CGL2005-00371/BOS, 
CGL2008-00378/BOS and CGL201 1-23466 to M.R. We thank P. Aguilera, N. 
Bonada, M. Charni, K. Gritzalis, A. H. Harrath, C. Hernando, L. Leria, C. Maritur, 
E. Mateos, M. Pala, S. Pons, I. Ribera, R. Sluys, E. Sola, J. Sola, G. A. Stocchino, 
M. Vila-Farre and M. Yacoubi for their collaboration on sample collection. G. 
A. Stocchino also for sharing information on the likely presence of D. sicula 
in Itlay. J. Rozas for his advice on the analytical methodologies. We also thank 
the Serveis Cientificotecnics (Universitat de Barcelona). Jim Hesson of Acade- 
micEnglishSolutions.com proofread the English. 

Received: 4 October 2013 Accepted: 3 December 2013 
Published: 11 December 2013 

References 

1. Lepori NG: Descrizione di Dugesia sicula, nuova sp. di Triclade d'acqua 
dolce dei dintorni di Catania. Arch Zool Ital 1948, 33:461-472. 

2. Benazzi M: Problemi di zoogeografia tirrenica studiati nelle Planarie. 

Atti Soc Tosc Sci Nat, Mem, (Serie B) 1 950, 58:21-28. 

3. Gourbault N: The karyotypes of Dugesia species from Spain (Turbellaria, 
Tricladida). Hydrobiologia 1981, 84(1 ):45-52. 



Lazaro and Riutort BMC Evolutionary Biology 2013, 13:268 
http://www.biomedcentral.com/1471 -21 48/1 3/268 



4. Pala M, Vacca RA, Casu S, Stocchino S: The freshwater planarian Dugesia 
sicula Lepori from Sardinia (Platyhelminthes, Tricladida). Hydrobiologia 

1995,310:151-156. 

5. Charni M, Ammar AB, Jaafoura MH, Zghal F, Tekaya S, Suppl: Spermatogenesis 
and spermatozoon ultrastructure in Dugesio sicula Lepori, 1948 
(Platyhelminthes, Tricladida, Paludicola). Belg J Zool 201 0, 1 40:1 1 8-1 25. 

6. Charni M, Ammar AB, Jaafoura MH, Zghal F, Tekaya S, Suppl: Ultrastructure 
of germaria and vitellaria in Dugesia sicula Lepori, 1948 (Platyhelminthes, 
Tricladida, Paludicola). Belg J Zool 2010, 140:1 11-118. 

7. Harrath AH, Sluys R, Merzoug D, Yacoubi-Khebiza M, Alwasal SH, Riutort M: 
Freshwater planarians (Platyhelminthes, Tricladida) from the Paleartic 
section of the African continent: new records, with the description of a 
new species. Zootoxo 2012, 3182:1-15. 

8. Lazaro EM, Sluys R, Pala M, Stocchino GA, Baguha J, Riutort M: Molecular 
barcoding and phylogeography of sexual and asexual freshwater 
planarians of the genus Dugesia in the Western Mediterranean 
(Platyhelminthes, Tricladida, Dugesiidae). Mol Phylogenet Evol 2009, 
52(3):835-845. 

9. Ribas M, Pala M, Vacca RA, Riutort M, Baguha J: Taxonomical status of the 
western Mediterranean asexual populations of the Dugesia (D) 
gonocephala group. Morphological, karyological and biochemical data in 
free-living and symbiotic Platyhelminthes. In Free-living and symbiotic 
Platyhelminthes. Fortschritte der Zoologie/Progress in Zoology Edited by Ax P, 
Ehlers U, Sopott-Ehlers B. Stuttgart: Gustav Fischer Verlag; 1988:129-137. 

10. Charni M, Harrath AH, Sluys R, Tekaya S, Zghal F: The freshwater planarian 
Dugesia sicula Lepori, 1948 (Platyhelminthes, Tricladida) in Tunisia: 
ecology, karyology, and morphology. Hydrobiologia 2004, 51 7(1 ):1 61 -170. 

1 1 . Stocchino GA, Manconi R: Overview of life cycles in model species of the 
genus Dugesia (Platyhelminthes: Tricladida). Ital J Zool 2013, 80(3)319-328. 

1 2. De Vries EJ: On the identity and occurrence of Dugesia sicula and D. 
biblica (Platyhelminthes, Tricladida, Paludicola) in the Mediterranean 
region. Fortschr Zool 1 988, 36:405-41 1 . 

13. Brown RM, Linkem CW, Siler CD, Sukumaran J, Esselstyn JA, Diesmos AC, 
Iskandar DT, Bickford D, Evans BJ, McGuire JA, Grismer L, Supriatna J, 
Andayani N: Phylogeography and historical demography of Polypedates 
leucomystax in the islands of Indonesia and the Philippines: Evidence 
for recent human-mediated range expansion? Mol Phylogenet Evol 2010, 
57(2):598-619. 

14. Lazaro EM, Harrath AH, Stocchino GA, Pala M, Baguna J, Riutort M: 
Schmidtea mediterranea phylogeography: an old species surviving on a 
few Mediterranean islands? BMC Evol Biol 201 1, 1 1(1):274. 

15. Zhang Y, Zhang D, Barrett SCH: Genetic uniformity characterizes the 
invasive spread of water hyacinth (Eichhornia crassipes), a clonal aquatic 
plant. Mol Ecol 201 0, 1 9(9):1 774-1 786. 

16. Ball IR, Fernado CH: Freshwater Triclads (Platyhelminthes, Turbellaria) and 
continental drift. Nature 1969, 221:1 143-1 144. 

1 7. Ball IR: Nature and formulation of biogeographical hypotheses. Syst Zool 
1975, 24:407-430. 

18. Gourbault N: Expansion de Dugesia tigrina (Girard), planaire americaine 
introduite en Europe. Ann Limnol 1969, 5:3-7. 

19. Ball IR: Dugesia lugubris (Tricladida: Paludicola), a European immigrant 
into North American fresh waters. J Fish Res Board Can 1969, 26:221-228. 

20. Sola E, Sluys R, Gritzalis K, Riutort M: Fluvial basin history in the 
northeastern Mediterranean region underlies dispersal and speciation 
patterns in the genus Dugesia (Platyhelminthes, Tricladida, Dugesiidae). 
Mol Phylogenet Evol 2013, 66(3):877-888. 

21 . Reddien PW, Sanchez-Alvarado A: Fundamentals of planarian regeneration. 
Annu Rev Cell Dev Biol 2004, 20:725-757. 

22. Salo E, Abril JF, Adell T, Cebria F, Eckelt K, Fernandez-Taboada E, 
Handberg-Thorsager M, Iglesias M, Molina MD, Rodnguez-Esteban G: Planar- 
ian regeneration: achievements and future directions after 20 years of 
research. Int J Dev Biol 2009, 53:1 31 7-1 327. 

23. Kmiec B, Woloszynska M, Janska H: Heteroplasmy as a common state of 
mitochondrial genetic information in plants and animals. Curr Genet 
2006, 50(3):149-159. 

24. Calloway CD, Reynolds RL, Herrin GL, Anderson WW: The frequency of 
heteroplasmy in the HVII region of mtDNA differs across tissue types 
and increases with age. Am J Hum Genet 2000, 66(4): 1 384-1 397. 

25. Bessho Y, Ohama T, Osawa S: Planarian mitochondria I. Heterogeneity of 
cytochrome c oxidase subunit I gene sequences in the freshwater 
planarian, Dugesia japonica. J Mol Evol 1992, 34:324-330. 



Page 13 of 14 



26. Hall TA: BioEdit: a user-friendly biological sequence alignment editor and 
analysis program for Windows 95/98/NT. Nucl Acids Symp Ser 1 999, 41 :95-98. 

27. Alvarez-Presas M, Carbayo F, Rozas J, Riutort M: Land planarians 
(Platyhelminthes) as a model organism for fine-scale phylogeographic 
studies: understanding patterns of biodiversity in the Brazilian Atlantic 
Forest hotspot. J Evol Biol 201 1, 24(4):887-896. 

28. Savolainen P, Arvestad L, Lundeberg J: mtDNA tandem repeats in 
domestic dogs and wolves: mutation mechanism studied by analysis of 
the sequence of imperfect repeats. Mol Biol Evol 2000, 1 7(4):474-488. 

29. Librado P, Rozas J: DnaSP v5: a software for comprehensive analysis of 
DNA polymorphism data. Bioinformatics 2009, 25(1 1 ):1 45 1 -1 452. 

30. Clement M, Posada D, Crandall KA: TCS: a computer program to estimate 
gene genealogies. Mol Ecol 2000, 9:1657-1659. 

31. Phillips SJ, Dudik M, Schapire RE: A maximum entropy approach to 
species distribution modeling. In ICML '04 Proceedings of the twenty-first 
international conference on Machine learning. ACM New York, NY, USA 
©2004; 2004:655-662. doi:1 0.1 1 45/1 01 5330.1 01 541 2. ISBN 1-58113-838-5. 

32. Phillips SJ, Anderson RP, Schapire RE: Maximum entropy modeling of 
species geographic distributions. Ecol Model 2006, 190(3-4):23 1-259. 

33. Elith J, Phillips SJ, Hastie T, Dudik M, Chee YE, Yates CJ: A statistical 
explanation of MaxEnt for ecologists. Divers Distrib 201 1, 17(1):43-57. 

34. Mourn T, Bakke I: Mitochondrial control region structure and single site 
heteroplasmy in the razorbill (Alca torda; Aves). Curr Genet 2001, 39(3): 198-203. 

35. Garcia-Diaz A, Oya R, Sanchez A, Luque F: Effect of prolonged vegetative 
reproduction of olive tree cultivars {Olea europaea L.) in mitochondrial 
homoplasmy and heteroplasmy. Genome 2003, 46(3):377-381. 

36. Williams ST, Knowlton N: Mitochondrial pseudogenes are pervasive and 
often insidious in the snapping shrimp genus Alpheus. Mol Biol Evol 2001, 
18(8):1 484-1 493. 

37. Richly E, Leister D: NUMTs in sequenced eukaryotic genomes. Mol Biol 
Evol 2004, 21(6):1081-1084. 

38. Song H, Buhay JE, Whiting MF, Crandall KA: Many species in one: DNA 
barcoding overestimates the number of species when nuclear 
mitochondrial pseudogenes are coamplified. Proc Natl Acad Sci 2008, 
105(36):! 3486-1 3491. 

39. Brabec J, Scholz T, Kralova-Hromadova I, Bazsalovicsova E, Olson PD: 
Substitution saturation and nuclear paralogs of commonly employed 
phylogenetic markers in the Caryophyllidea, an unusual group of 
non-segmented tapeworms (Platyhelminthes). Int J Parasitol 2012, 
42(3):259-267. 

40. Bogenhagen D, Clayton DA: The number of mitochondrial 
deoxyribonucleic acid genomes in mouse L and human HeLa cells: 
quantitative isolation of mitochondrial deoxyribonucleic acid. J Biol Chem 
1974, 249(24):7991-7995. 

41. Shadel GS, Clayton DA: Mitochondrial DNA maintenance in vertebrates. 
Annu Rev Biochem 1 997, 66:409-435. 

42. Walton C, Butlin RK, Monk KA: A phylogeny for grasshoppers of the genus 
Chitaura (Orthoptera: Acrididae) from Sulawesi, Indonesia, based on 
mitochondrial DNA sequence data. Biol J Linn Soc 1997, 62(3)365-382. 

43. Sachadyn P, Zhang X, Clark LD, Naviaux RK, Heber-Katz E: Naturally 
occurring mitochondrial DNA heteroplasmy in the MRL mouse. 
Mitochondrion 2008, 8(5-6)358-366. 

44. Janko K, Drozd P, Eisner J: Do clones degenerate over time? Explaining 
the genetic variability of asexuals through population genetic models. 
Biol Direct 201 1, 6(1):17. 

45. De Vries EJ, Ball IR: On Dugesia gonocephala from Western Europe. Contrib 
Zool 1980, 50(2)342-350. 

46. Carbayo F, Alvarez-Presas M, Olivares CT, Marques FPL, Froehlich EM, Riutort 
M: Molecular phylogeny of Geoplaninae (Platyhelminthes) challenges 
current classification: proposal of taxonomic actions. Zoologica Scripta 
2013, 42:508-528. 

47. Birky- CW Jr: Heterozygosity, heteromorphy, and phylogenetic trees in 
asexual eukaryotes. Genetics 1996, 144(1)427-437. 

48. Mark Welch DB, Meselson M: Evidence for the evolution of bdelloid 
rotifers without sexual reproduction or genetic exchange. Science 2000, 
288(5469):1 21 1-1215. 

49. Stocchino GA, Corso G, Manconi R, Pala M: African planarians: Dugesia 
aethiopica sp. n. (Platyhelminthes, Tricladida) from Lake Tana 

(NW Ethiopia). Ital J Zool 2002, 69:45-51. 

50. Rosenbaum G, Lister GS: Formation of arcuate orogenic belts in the 
western Mediterranean region. In Orogenic Curvature: Integrating 



Lazaro and Riutort BMC Evolutionary Biology 2013, 13:268 
http://www.biomedcentral.com/1471 -21 48/1 3/268 



Page 14 of 14 



Paleomagnetic and Structural Analyses. Edited by Sussman AJ, Weil AB. 
Boulder, Colorado: Geological Society of America; 2004:41-56. 

51 . Schettino A, Turco E: Plate kinematics of the Western Mediterranean 
region during the Oligocene and Early Miocene. Geophys J Int 2006, 
166(3):1 398-1 423. 

52. Calow P, Beveridge M, Sibly R: Heads and tails; adaptational aspects of 
asexual reproduction in freshwater triclads. Am Zool 1979, 19:715-728. 

53. Vila-Farre M: Planarias ibericas: taxonomi'a y biogeografia. In PhD thesis. 
2010: Universitat de Barcelona, Departament de Genetica; 2010. 

54. Cubeddu T, Alberti A, Manconi R, Pala M: Studio di una popolazione 
scissipara di Dugesio siculo Lepori (Platyhelminthes, Tricladida) dell'isola 
di S Antioco. Biogeographia 1995, 17:193-197. 

55. De Vries EJ: Further contribution to the taxonomy and biogeography of 
the subgenus Dugesia (Platyhelminthes: Tricladida: Paludicola) in the 
Mediterranean region and the Middle East. Isr J Zool 1988, 35:109-136. 

56. Reynoldson TB: Observations on the freshwater triclads of North Wales. 
Ann Mag nat Hist SI 2 1 956, 9:61 2-622. 

57. Wright JF: Colonization of rivers and canals in Great Britain by Dugesia 
tigrina (Girard) (Platyhelminthes: Tricladida). Freshwat Biol 1987, 
17(1):69-78. 

58. Young JO, Reynoldson TB: Continuing dispersal of freshwater triclads 
(Platyhelminthes; Turbellaria) in Britain with particular reference to lakes. 

FreshwBiol 1999, 42:247-262. 

59. Batistoni R, Filippi C, Salvetti A, Cardelli M, Deri P: Repeated DNA elements 
in planarians of the Dugesio gonocephalo group (Platyhelminthes, 
Tricladida). Hydrobiologia 1998, 383(1 -3):1 39-1 46. 

60. Batistoni R, Rossi L, Salvetti A, Deri P: A molecular cytogenetic comparison 
of planarians from the 'Dugesia gonocephala group' (Platyhelminthes, 
Tricladida). Ital J Zool 1999, 66:239-244. 

61. Benazzi M, Deri P: Taxonomic perspectives concerning fissiparous 
populations of the planarian Dugesio gonocephalo s.l. inferred from 
ex-fissiparous specimens. Atti Accademia Nazionale del Lincei 1988, 
19:45-64. 

62. Deri P, Colognato R, Rossi L, Salvetti A, Batistoni R: A karyological study on 
populations of Dugesio gonocephalo s.l. (Turbellaria, Tricladida). Ital J Zool 
1999, 66(3):245-253. 

63. Bromley HJ: Morpho-karyological types of Dugesio (Turbellaria, Tricladida) 
in Israel and their distribution patterns. Zool Scr 1974, 3(5-6):239-242. 

64. Bromley HJ: Observations on fission, maturation and reproduction in 
Dugesio biblico (Turbellaria, Tricladida) from Israel. Zoologica Scripta 1977, 
6:197-201. 

65. Reynoldson TB: The quantitative ecology of lakedwelling triclads in 
northern Britain. Oikos 1958, 9:94-138. 

66. Reynoldson TB, Jaques RM: Changes in the triclad (Turbellaria, 
Platyhelminthes) fauna of some North Wales lakes during the 20 years 
up to 1 973. Nature in Wales 1 976, 1 5:73-80. 



doi:1 0.1 186/1471-2148-13-268 

Cite this article as: Lazaro and Riutort: Dugesia sicula (Platyhelminthes, 
Tricladida): the colonizing success of an asexual Planarian. BMC 

Evolutionary Biology 201 3 1 3:268. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at (^\ RiftMM i rpntral 

www.biomedcentral.com/submit \^ ™omea centra I 



